arXiv:1508.05666v2 [hep-th] 30 Oct 2015 


Path-Integral Approach to the Scale Anomaly at Finite Temperatnre 


Chris L. Lin and Carlos R. Ordonez 
Department of Physics, University of Houston, Houston, TX 77204-500^ 

(Dated: November 2, 2015) 

We derive the relativistic thermodynamic scale equation using imaginary-time path integrals, with 
complex scalar field theory taken as a concrete example. We use Fujikawa’s method to derive the 
scaling anomaly for this system using a matrix regulator. We make a general scaling argument 
to show how for anomalous systems, the /? function of the vacuum theory can be derived from 
measurement of macroscopic thermodynamic parameters. 
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I. INTRODUCTION 

In a series of seminal papers by Callan, Coleman, and Jackiw 0,1, it was noted that in general the trace of the 
Belifante stress-energy tensor 0^) for any renormalizable theory could be improved, so that classically for scale-invariant 
systems (systems invariant under the conformal group). 


01 ) = 0 . ( 1 ) 

This improved tensor has a number of desirable properties over the canonical tensor (the one derived from Noether’s 
theorem) such as having finite matrix elements in the quantum theory, and that the energy for bound states can 
be naturally expressed as the trace of this tensor. Shortly after these observations, it was noted that the same 
improvement program could be applied in the non-relativisitic case Q, so that for classical scale-invariant systems 
(systems invariant under the Schrbdinger group): 


3 

20°°-^0"=O, (2) 

i=l 

where the 2 results from the fact that in non-relativisitic theories time must scale as twice the power of space. ^ 

Eqs. (HI) and ([2]) fail to consider the trace anomaly. In general, the trace of the stress-energy tensor taken between 
bound states gives the energy of the bound state: 


Ef,= I dU(0p, (3) 

which derives from the fact that the time average of the field virial is zero for bound states 0. With slight 
modification Eq. m holds in the non-relativistic case too (see 0 for a specific example). However, it is well known 
that even though 0^=0 for a classically scale-invariant system, which would imply bound states can only have zero 
energy,^ the quantization procedure can destroy this relationship. When this happens this is called a scale anomaly, 
and is the mechanism that allows the bound state energy to differ from zero. 

As an example, in QCD with massless quarks (or no quarks at all), the Lagrangian is classically scale invariant so that 
= 0. However, through the renormalization process, a scale appears as Aqcd- In general this makes (0()) = A, 
where A is the anomaly. The stress-energy tensor can then be further improved: 
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^ The Schrodinger equation has only one derivative of time, and two of space, so for scale invariance time must scale as twice the power 
of space. 

^ This is also obvious from the fact that there are no scales to even form Ef,. 



so that is no longer traceless. Then 


Eb = J dV (T™) = J dV (9°°) + (5) 

which implies that A accounts for 1/4 of the energy of the hadron. This can explicitly be seen in the bag model where 
confinement of the quarks and gluons is the result of a cosmological constant term in the Lagrangian which contributes 
a positive energy and negative pressure to 9^^", which confines the system. Then from the tracelessness of 9^’', 
A = so that confinement accounts for 1/4 of the hadron energy [^. 

In this paper, we are interested in the thermal analogues of Eqns. dm and Both of these quantities are very 
important in their respective areas of physics. In the nonrelativistic sector, for an ultracold dilute gas, (HI) would read: 

j-2 

2£ — 2,P=—- — X{{tp\x)ijj{x)Y) . ( 6 ) 

The RHS is known as the Tan contact, and is extremely important in atomic physics. In terms of it, Tan derived a 
set of universal relations that govern many relationships between the thermodynamics variables of the system 
and the behavior of the large momentum tails of correlation functions. These relationships hold even in the strongly 
interacting regime where perturbation theory becomes inadequate [B]. A field theoretic explanation of Tan’s result 
was later developed in terms of the operator product expansion la¬ 
in QCD, the analog would be 0: 
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£-W = Yy, . (7) 

In the low temperature regime where the coupling g is strong, the trace anomaly of the RHS is calculated by calculating 
the LHS of Eq. ([^ using a lattice action. The goal is to calculate the QCD equation of state P = P{T,fi, V) rather 
than the anomaly itself. However, for technical reasons [iB]) ^ ~ is important as an intermediate step in lattice 
QCD for calculating P{T,fi, V), where it is given by: 


A = £- 


3P = - 


T dlnZ 
Edina’ 


and plugging into Eq. 0 gives after using thermodynamic identities: 


( 8 ) 


d / P 
dlnT 

which can then be integrated to get P{T, g,^V). a is the lattice spacing and Z is the partition function with lattice 
action. 

In this paper, following the approach initiated in [isl - fTsjl for non-relativisitic systems, we provide a continuum/non- 
lattice path-integral approach to deriving the thermodynamic trace equation £ — SP, where anomalies naturally appear 
as a result of a change of variables of the path-integral measure, the thermal analog of Fujikawa’s method. This is in 
contrast to an operator approach, where one takes the thermal quantum statistical expectation values of both sides 
of Eqns. like o and and identifying = £ and = Ph, where Ph is the hydrodynamic pressure [lB| . 

Within this path-integral approach, no reference needs to be made about improvement of the stress-energy tensor, or 
the validity of equating the hydrodynamic pressure Ph with the thermodynamic pressure P derived from the grand 
partition function, which is nontrivial, especially in the presence of anomalies |17l. Il8l|. For concreteness, we will take 
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as our system a complex scalar field theory, but the results can be extended for other systems. The Lagrangian is 
given by 


£ = — rn^cji^cj) — 


( 10 ) 


and has a 17(1) symmetry 


leading to a conserved charge: 


Under scale transformation: 


(f) 


jo = 

Q = i [ d^x 4>^do4>- 


(j)'{x') = e 
(j)'^{x') = e~^(j)^{x). 


( 11 ) 


( 12 ) 


(13) 


II. THERMODYNAMIC DILATION EQUATION 

For a homogeneous system the grand potential Q = Q(/3,/i, U) in the large volume limit equals —PV, so that the 
partition function is Z = , and can be expressed via a path integral: 




(14) 


with^ 


Se= [ [ d^xdr (d^(j}*d^(l)+, 

Jo Jv \ 4 y (15) 

jo = -(j)*dr(j). 

Now consider an infinitesimal “relativistic thermodynamic scaling” 


f3' = ePI3 = p + pf3 = l3 + 6(3, 

L( = e^Li = Li + pLi = Li + SLi, (16) 

p' = p. 

where Li is the length of the box in the i direction and p is a dimensionless infinitesimal parameter. 
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Due to the dependence of j-'o on conjugate momenta, when integrating out conjugate momenta to pass into the Lagrangian formulation 
of the path integral, Ce acquires an additional term: see [T^l . 
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In the large volume limit it is assumed that P{f3, V) = P{I3, so under the transformation of Eq. (fT6l) : 


5{I3PV) = {6f3)PV + I3{SP)V + PP{SV) 

= p(^ppV + p(^^P^V + pPi3V)y 

Now using the identity fiV = —PV — E + fiQ, we get 

SipPV)=pi-pE + pPi3V)+(3pQ), (18) 

and therefore 


5 = d il3PV) = p/S {-E + 3PV + pQ) . (19) 

Eq. (jiyp represents the effect of the scaling in Eq. (fTCl) on the LHS of Eq. (1141) . Now we analyze the effect of this 
scaling to the RHS of Eq. da, the path-integral part, from which anomalies originate, and eventually equate the 
two expressions. 

The scaling in Eq. da represents a dilation of the system: 


e P(j){x), (20) 

The dilated system has 


<j>'{x') = 

r{x) = 


where 


^P'P'V _ 




( 21 ) 


^E = y d^x'dr' + {w? - -h ^(<()'* 


rP' r /■e'’/3 /• / <_). 

p / d^x'dT'jQ = p / / d^x'dr' \ 

Jo Jv Jo JePV 


( 22 ) 


To compare to the undilated system, we “pull back” to unprimed variables by substituting Eq. (1201) into Eq. dUD 
and Eq. (1^ . Eq. (1^ becomes: 


S'e = £ d^x'dr' + irn^ - p^)c^'*4>' + ^( 0 '* 

= Se + 2p f j d^xdr {mJ — p'^)(j)*(j). 

Jo Jv 


A 


+ :-(e-2p<),*0)2 


(23) 


^ This can be shown via cluster decomposition: e.g., see M- 
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Similarly; 


fP' f /■e'’/3 

H / (fx'dr' jo = ^J. 

Jo Jv' Jo 


d^x'dr' -(h'*d: 


'epy 

73. 


(24) 


=/r / / d xdr jo + pfi / / d xdrjo. 

Jo Jv Jo Jv 


Plugging in these expressions into Eq. (EU; 


g/3'P'y' _ J d?xdT {Tn?-p'^)!j>’’ 4 ,+pp J^d^xdTjo^ ^25) 

where J is the Jacobian of the transformation (<()', (j)'*) —>■ (c^, (j*). Expressing J = 1 — pA and using Eq. ([T9]): ® 


J [ef^PV] = p/3 {-E + 3PE + pQ) = p d^xdr (' 

The chemical potential terms drop out on both sides® and we get: 


- p"‘)(t)^4> ) ^y J e' 


ppv 


(26) 


where 


J = 


dfj'dc/)'' 


£-3P = 2m‘^ {(jU) + A 


— gTrlog(72(<5^(a:-y)+p(-l-a:^a^)<5^(a:-y))) 


d(j)d(j)* 

_ p f d'^xti:[{-l-Xij,d^)S‘^(x-y)l2]\ 


SO that 


= l + p J d'^xtr [{-l-x^df,)S'^{x-y)l 2 ]\^^^ 


A = tr [(1 + Xf,df,)6'^{x - y)l 2 ] 


(27) 


(28) 


(29) 


I 2 is the two dimensional identity matrix which results from having two fields, (j) and (j)* 7 A = is the anomaly, a 
divergent quantity that requires regularization. 


III. FUJIKAWA CALCULATION 

In Euclidean space, Le = d^(j)'d^(j)+'>A4>A+j{A4'Y- A saddle point expansion about a constant classical background 
(j) produces the quadratic piece C 2 ' 


f, 1/t \ (A rr? E ^ 


\ fv 


- ^ (ut A (-9^ + ^ I# V ^ 


( 30 ) 




® i I[d/>][d/>*]F(<j},<j}*)e~^‘^P'"^o Jyd^xdTjo^ 

® Using the identity Q = and Eq. (I14II . Q = = {jo) + 2fi 

^ Note that Tr in Eq. (I28I I refers to both discrete (2 x 2) and continuous variables, whereas tr in Eq. [12911 refers to only (2 x 2). 
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where C = + X(j)*(j), Tj is the fluctuating field around (j), and M is a Hermitian matrix. Following Fujikawa [^ . 

we use M, the bilinear matrix, as the Hermitian matrix that goes in our regulator®. Choose a regulator of the form 
R = R[^) with the property that R{0) = l 2 - The expression to be regulated is: 


where 9 = 1 + Xudn,, so that 


This expression equals: 


^ = tr 


65{x — y) 0 
0 95{x — y) 


(31) 


Ar = tr 


■^ ( ^ ) SS{x-y)l2 


(32) 


I -o'Ac 


= A4 
= A4 


d^k 

W 

r d'^k 


2A2 

k^+C 


X(f)<p 


i^c I 


2h? 


A2 


C X(p(p 
tr 7? I ',7^ 2A' 

(2^)4 


I \ 

) (1 - iAxfj.kfj,) 
A^/ 


( 2 *)“ V ^ (,.2 + 4 


A2 


(33) 


where the term is odd so vanishes over the integral when multiplied by the even function R(—k) = R{k) = f{k^). 
Next we define: 


D = k^h, 

I ( C (34) 

« - A2 c ) ’ 

so that the equation can be written succinctly: 

Ar=A^ J ^tTR{D + B). (35) 

We then Taylor expand about D (note that [D,B] = 0 so the Taylor expansion is valid): 

Ar = A^I ^tr (^R{D) + R'{D)B + \r"{D)B^ + ..}j . (36) 

The first term is the same as in the non-interacting case, which is taken to be anomaly free [2l|, so we neglect it. The 
second term can be absorbed by a mass counterterm. Terms higher order than the third term fall faster than ^ so 
the A^ prefactor in Eq. (1361) cannot keep them from going to zero. Only the 3rd term is independent of the cutoff. 
Therefore: 


e.g., for the chiral anomaly with C = the matrix ilj) is to be used as the argument of the regulator. M, the quadratic piece of 

the quantum action, naturally captures the 1-loop effects of interactions which are responsible for anomalies. 
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where the solid angle Vl = 2 tt'^ was used. Now 


( 37 ) 


r.2^±(c^ + AC# \ ^ B2\ 

\ XC(l)*(j>* + A!(£1^ J A2 Bi) 

and since R[D) is diagonal, we can define: 


(38) 


R{D) = f{e)h. 

Note that the derivative in Eq. (|37| is w.r.t. . Therefore: 


, . 4 f k'^dk^ 1 ! n , ^ 

Ar = A'^ I ^^-tr{R"{D)B^ 


= B, 


16# 2 

k^dk"^ 


IGtt^ 


rm, 


where we have safely taken A —^ oo. Integrating by parts: 


(39) 


(40) 


where we require 


Ar = 


IGtt^ 

Bi 

IGtt^ 

Bi 

IGtt^ ’ 


[fc#'(#)] 

[#/'(#)] 



/(O) = 1 

/(oo) = 0 

[#/'(#)] 1“ = 0, 

which are the same conditions on the regulator for the chiral case [^ . 
Plugging in Bi from Eq. (l38l) into Eq. (ITO . we get: 


(41) 


(42) 


Ar 


(j2 _j_ 

IOtt^ 


5A^(0*0)^ rn^ Am^ (/>*/>) 
647r^ IGtt^ Stt^ 


(43) 


The second term is independent of the coupling, and since the free theory is taken to be non-anomalous, we can 
subtract it. The third term can be absorbed into the mass term of Eq. (Eli), leaving only the 1st term as the anomaly 
[ 2 ^. Therefore 


Note that the anomaly Ar occurs inside the path integral, and ^/[##*]/(# (/)*)e“‘^®+ - = {f{4>,4>^)), so that in 
Eq. (ITil) there are expectation values. This replacement is valid up to 1-loop [^ . 






















IV. DIMENSIONAL ANALYSIS FOR RELATIVISTIC SYSTEMS 


In relativistic theories we set fi = c = ks = ^- The units for all quantities can then be written as fi'c^kgL^, where L 
is a variable in the problem with units of length. Suppose the system has microscopic parameters gk, which can be 
coupling constants or dimensionally transmuted quantities. We define [g^] = i as the power of L when g^. is written 
in units of Tic^k^L^. So for example [m\ = [E] = —1. The grand potential 0 = VL{j3, gi^V^gi) has [O] = —1 and can 
be written as: 


n{P,z,,V,gi) = (45) 

where /(^i,jg g, dimensionless function of dimensionless variables, Zi is the fugacity corresponding to /i^ 
{zi = and D is the number of spatial dimensions.® VL has this form because /3 and /ii don’t depend on the 

absolute size of the system (they are intensive variables). If one doubles the system keeping /3 and gi constant, then 
being an extensive quantity, should double. So il must be proportional to V.^® To make up for the remaining 
dimension ([fl] = —1), we are free to pull out one of the dimensionful arguments of f2, and the rest of the arguments 
must be ratios with the argument we pulled out. We will pull out /3. This is equivalent to choosing our scale as /3 
and measuring all other quantities in units of (3. 

Take the derivative of Eq. (I45p w.r.t. to f3 at constant fugacity Zi and volume V, and multiply times /3: 




= {-1- D)il + VP 


Zi,V 


-1-D , 


df{z,,g,/3-^s.]) 


-1-D f 


= {-1- D)n + vp 


= D)n-Y^[gk]gi 


E 


dp 


-[ 9 k]gk df{zi,giP [94) 


P 


dgk 


an 

dgk' 


(46) 


Now, we use the thermodynamic identity E = 


Zi,V 


— n + /3 


Zi,V 


E-DPV= 


- DPV 


Zi,V / 


an 


= + (-1 - Z4) j - DPV 

= -{^P+{-\-D)P- V - DPV 


f - DP = y^[gfc]gfc 


dP 

dgk' 


(47) 


where the derivatives are at constant /5, and V. 


^ For example, if the coupling gi has dimensions of length, the corresponding dimensionless variable is gip~^ = giT which is dimensionless. 
If the coupling g 2 as dimensions of energy, = 920 = 

D, = —PV, so Eq. (I45II is consistent with the statement that P(/3,/j,, V") = P(0,g,). 
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FIG. 1: Diagrams contributing to the P function for complex scalar field theory. 1 and 2 refer to incoming particles, 3 and 4 
to outgoing particles. 


V. P FUNCTION 


For a system that develops a microscopic scale M through dimensional transmutation via renormalization of the 
coupling constant: 


dM dX dM dX 



(48) 


since If = pulls 

thermal average. 


down the interaction term in the path integral, creating a 


Comparison of Eq. (0, Eq. (I?7l) . Eq. (H^ . and Eq. (H51) gives: 


as 


m 


5A" 
IGtt^ ’ 


(49) 


would give Eq. pg]). 

The p function of Eq. (|4^ can be gotten from setting e = 0 for the charge e in the calculation for the four-scalar 
vertex in scalar electrodynamics [2^. A diagrammatic calculation requires the identification of 3 diagrams (see Fig[T]). 
Diagram (a) contains a symmetry factor of 1/2 due to the swapping of internal propagators. Modulo the symmetry 
factor, each diagram contributes the same amount to the P function, giving 1/2 -|- 1 -|- 1 = 5(1/2), or the first 
diagram’s contribution multiplied by 5. The matrix M used for regularization automatically mixes the interactions, 
giving the factor of 5. Using the definition of the beta function and setting the renormalization scale 

M = T, one can solve the differential equation for the coupling X{T) = v , where A is the Landau pole. As 

f —>• 0 the coupling is small and the system behaves like a gas of noninteracting bosons, while as T —A the coupling 
blows up and perturbation theory fails. 


VI. CONCLUSIONS 

In this paper we have extended to relativistic systems the path-integral approach to the study of quantum anomalies 
for many-body systems initiated in [IM3. A notable difference is that in the relativisitic case we have a very wide 
class of regulators characterized by the function f{k^) of Eq. (IMl) . which other than satisfying Eq. (1421) . is of a very 
general nature. An interesting result of this paper is the extraction of the leading order result for the beta function 
for complex helds, Eq. (|4^ . obtained here by comparing Eqs. (Uni), dUl), dill) and (l48l) , without resorting to graphical 
methods [^ . This result gives further support to the importance of Fujikawa’s approach in the description of 
quantum anomalies for systems at finite temperature and density. We are currently pursuing further studies and 
extensions of this method, as well as applications to other systems with classical scale symmetry. 
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